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Abstract 

We present expressions for the coefficients which arise in asymptotic expansions of multiple integrals 
of Laplace type (the first term of which is known as Laplace's approximation) in terms of asymptotic 
series of the functions in the integrand. Our most general result assumes no smoothness of the functions 
of the integrand, but the expressions we obtain contain integrals which may be difficult to evaluate in 
practice. We then make additional assumptions which are sufficient to simplify these integrals, in some 
cases obtaining explicit formulae for the coefficients in the asymptotic expansions. 
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1 Introduction. 

Consider the integral 

b 

e- kf(x) g(x)dx 

where / and g are sufficiently smooth functions and / attains a nondegenerate, unique minimum in [a, 6] at a 
point xq G (a, b). In Lap95o Laplace observed that as k increases, the integral localizes near the minimum 
of /, and in particular showed that as k —> oo, 

In this article we will be concerned with the asymptotics of integrals of the form 

JR 

known as Laplace-type integrals, where k is a real parameter which goes to infinity, / and g are real- valued 
functions defined on some region R C M d , and we assume that / attains a unique minimum in the interior 
of R. For convenience, say this minimum value of / is at € R. 

Laplace's idea remains valid in d dimensions; namely, for k sufficiently large, the integral becomes localized 
near the minimum of / and one may compute what is now known as Laplace's approximation: 



R 



;Wx ^r *<s> (L2) 



k J y/dctHf(0) 



where H f(0) denotes the Hessian of / evaluated at 0. To prove this requires two main ideas: first, one shows 
that the integral localizes to a neighborhood of 0. Second, one applies the Morse Lemma to transform /, and 
hence the integral, to a standard form from which one can deduce the desired approximation. We refer the 
interested reader to [BH75| . |dB81) . and jWonOl] for some modern treatments of Laplace's approximation 
and related ideas. 
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1 The approximation HI. Il l appears in Lap95 , Part 2, Chapter 1, Section 27. 



1 



When the function / attains a unique minimum somewhere on the boundary, the integral localizes near 
this boundary point as k — > oo, but the analysis and resulting asymptotics are slightly different. We do not 
consider this case in this article. 

There are several important and useful generalizations of Laplace's approximation. When / is purely 
imaginary, because of cancelations arising from oscillations, the integral localizes near critical points of / (not 
just minima), and the resulting Laplace- type approximations are called stationary phase approximations. If 
/ is defined on a region in C and is complex valued, the analysis takes on quite a different character, and 
Laplace- type approximations are known as the method of steepest descent. 

All of these Laplace-type approximations are very useful in applications. For example, Laplace and 
stationary phase approximations arise naturally in statistical mechanics and quantum field theory because 
the partition functions tend to be of Laplace type. In the quantum case, the natural parameter is i/h 
and one is interested in the semiclassical limit h — > 0, whereas in statistical mechanics (for example in the 
microcanonical ensemble) one is typically interested in the joint limit as the inverse temperature (which 
plays the role of k) and the dimension d both tend to infinity. 

Laplace-type approximations also play a role in pure mathematics; equivariant localization is an applica- 
tion of stationary phase to certain integrals of equivariant cohomology classes where the resulting asymptotic 
expansion is, somewhat magically, exact. In geometric quantization, Laplace's approximation has been used 
by the author and Brian C. Hall HK06. to analyze the non-unitarity of "quantization commutes with re- 
duction. "0 This list is by no means exhaustive or unbiased, and we urge the interested reader to consult the 
above-mentioned references for further details and applications. 

In this article, we will stay in the realm of real-valued functions and study a different extension of 
Laplace's approximation. In modern terms, Laplace's approximation gives the first term in an asymptotic 
expansion 

J e- k fgd d x^J2k- {j+d)/2 C^ 

This means, by definition, thal|f| 

/ e- k fgd d X = V fc-«+ d )/ 2 Ci + o(fc-( JV+d )/ 2 ). 

The main result of this paper is Theorem l2.1[ which gives expressions for Q under quite general hypothe- 
ses; in particular, we make no smoothness assumptions on / or g, assuming only that they admit asymptotic 
expansions as |x| — > 0. (The order of the expansions of / and g determines the order N of the asymptotic 
expansion of the integral.) The function g is even permitted a mild singularity at 0. The proof is essentially 
a modification of an existence proof of Fulks and Sather |FS61) using Theorem 12. 2[ which is a variant of a 
result of Frame |Fra57j regarding inversions of series. 

The expressions for the coefficient which appear in Theorem 12.11 though, involve certain integrals which 
may in general be difficult to evaluate. With more restrictive assumptions, we can apply Theorem !2.1l to go 
much further, obtaining Corollaries 13.21 and 13.31 and Theorem 1 1.1|) . These more restrictive assumptions are 
still sufficiently general to include most common applications (though certainly not all). In particular, we 
make two (partially) independent simplifications. 

The first simplification occurs if we assume that // |x|" is continuous at x = 0, where v > is the order 
of the zero of / at (c.f. Corollary I3.2j) . In this case, the integrals appearing in Theorem 12.11 simplify 
greatly, and given the appropriate data about / and g, one expects that the resulting integrals can be easily 
evaluated. 

The second simplification occurs if we assume instead that / has a nondegenerate minimum at 0. In 
this case, the Morse Lemma tells us that in a neighborhood of 0, the function / can be put into "standard" 
quadratic form. This extra structure allows us to make contact with the usual expressions of Laplace's 

2 "Quantization commutes with reduction" is a result of Guillemin and Sternberg GS82 which says that, at the level of 
vector spaces, the geometric quantization of the symplectic reduction of a Kahler G-manifold is isomorphic to the subspace of 
G-invariant vectors of the quantization of M itself. These spaces, though isomorphic, are not naturally unitary. 

3 A function h is little-o of k p , written h{k) = o(k p ), k — ¥ oo, if Wm^^^ k~ p h(k) = 0. A related notion which we will also 
use is big-O; h(k) = 0(k p ), k — > oo if there exists a constant C > such that for k sufficiently large, < Ck p . It follows 
that if h{k) = o(k p ), fc ^ oo, then h(k) = 0{k p ). Conversely, if h(k) = 0{k p+e ) for any e > 0, then h{k) = o(k p ). 
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approximation. Moreover, one can insure (by an appropriate linear change of coordinates) that //|x| 2 
is continuous at 0, and thus our previous simplification applies. The result is the following theorem (we 
emphasize that the main result of this paper is Theorem 12.11 we quote the following theorem because it 
involves only Taylor coefficients and combinatorial quantities, and so is in some sense our most "explicit" 
result, although it is simply a corollary of Theorem 12. II) . 

Theorem 1.1 Let R C M. d be a measurable set which contains as an interior point and suppose f E 
C N+2 (R) and g e C N (R). Suppose moreover that f has a unique, nondegenerate minimum value ofO at 0. 
Assume that for some fco > 0, the integral f R e~ k °f g d d x converges. Then there exists a linear transformation 
P : R d — > R d and an asymptotic expansion 

„ LW/2J 

/ e- kf gd d * = k- d / 2 V C2,/fc- J +0(£r (Ar+d)/2+1 ), fc^oo (1.3) 
where \_N/2\ denotes the largest integer less than N/2, and the coefficients are given by 

( 2 ^/ 2 * " ( _l)r 

x Y Y ■ ■ ■ Y even((3 + a x + ■ ■ ■ + a r ) 

\P\=2j— m !ni+n 2 -\ hn r =m\ |ai|=ni+2 \a r \=n r +2 

\ ni>l, i— l,...,r / 



x 



In/,!. 



in which H x f(0) is the Hessian of f at x = (with respect to x), y = Px are coordinates in which 
H y f(0) = 1]r™, and for a multi-index a — (a 1 , a 2 , . . . , a d ), Dyg denotes the mixed partial derivative of order 
a of g with respect to y, and 



even(a) 



if, for any m, a m is odd 

1 otherwise. 



Empty sums, 0! and (—1)!! are all understood to be 1. 

Recall that a critical point of / is said to be nondegenerate if the Hessian at that point is invertible. 
The coordinate change y = Px appearing in Theorem 11.11 can be computed explicitly by diagonalizing the 
Hessian of /. More precisely, since it is symmetric, we can write Hf(0) = Q~ 1 DQ for some orthogonal 
matrix Q and some diagonal matrix D. Since the minimum of / at is nondegenerate, the eigenvalues of 
Hf(0) are all positive. The matrix P of Theorem 1 1.1 1 is then P = \J~DQ. 

We give expanded expressions of the first few coefficients in the appendix. 

Our results are not the first of their kind. It has been known for some time that with sufficient hypotheses 
on / and g, such an expansion exists, and various methods for computing the coefficients have been given. 
In |BH75 ! , Bleistein and Handelsman assume that the minimum of / is nondegenerate to obtain a complete 
asymptotic expansion of the form (|1.3j) . The coefficients, though, are computed in terms of a Jacobian of 
a coordinate transformation (which essentially arises from the Morse Lemma), that in general cannot be 
computed explicitly. Nevertheless, since the result is eventually evaluated at x = 0, it is possible to proceed 
term by term, obtaining explicit formulas in terms of the derivatives of / and g. 

In a similar direction, in [Ski80] , assuming f{x) has a nondegenerate minimum, Skinner makes a coordi- 
nate change (which diagonalizes the Hessian of f(x) at the minimum) and computes £2- As he points out, 
his methods could likely be used to directly compute the higher coefficients in this case. The end result 
would be essentially Theorem 11.11 quoted above. 

In dimension one and assuming / = 'Y^jL2 a 3 x ' > > °2 7^ and g = Y^jLo^j xJ > de Bruijn |dB811 (4.4.9)] 
gives a complete asymptotic expansion of e~ k f g dx in terms of the coefficients of the series expansion 
of the product c/exp{— fcX^3 CLjX 3 }. This series expansion of the product can in principle be computed in 
terms of the dj and bj using Faa di Bruno's formula, though de Bruijn does not do this. 
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In the 1-dimcnsional case, Wojdylo gives a closed form expression for Q for all j Woj06a, Woj06b|. His 
expressions should correspond to de Bruijn's expansions in terms of Taylor series with the above mentioned 
simplifications using Faa di Bruno's formula. Our results are in fact a generalization of those of Wojdylo in 
the precise sense that when d=l, our formulas reduce to his (see Section . 

More recently, Denef and Sargos [DS92 , Kaminksi and Paris [KP9 8a, KP98b|, and Liakhovetski and Paris 
[LPOlj use the relationship between polynomials and Newton polygons to analyze the asymptotics when the 
exponent / is a polynomial, possibly with noninteger powers. In a related direction, Dostal and Gaveau 
DG89] also work with the Newton polygon to study stationary phase for a polynomial phase function. Our 
methods are not generally well adapted to the case of polynomials with noninteger powers since one must 
essentially make a Taylor expansion of the polynomial, at which point it is difficult to see the geometry of 
the Newton polygon. 

The rest of the paper is organized as follows. In Section [2] we prove our most general result, Theorem 12.11 
For the proof, we recall and prove a variant of a result of Frame |Fra57j on inversion of series, and summarize 
in Lemma 12.31 the results of Fulks and Sather |FS61] which we need. Section [5] concludes with a few brief 
remarks regarding some of the combinatorial quantities which arise in Theorem 12. II 

In Section I5TT1 we describe the simplifications that occur in Theorem 12.11 when we assume that // |x|" 
is continuous at for some (maximal) v > 0. Then, in addition to the continuity assumption at 0, we 
place additional smoothness hypotheses on / and g, and give a corollary to Theorem 12.11 which expresses 
the coefficients Q in terms of the Taylor coefficients of / and g. 

In Section [3.21 we suppose that the minimum of / is nondegenerate. First, we exploit the nondegeneracy 
of the minimum to show that by a linear change of coordinates, we can insure that // |y| 2 is continuous at 
0. We then use the results of Section [3~T1 to obtain Theorem ll.il Finally, we relate Theorem l2.1l to Laplace's 
approximation (ll.2[) . 

In Section I3~31 we show that in dimension one our results reduce to those of Wojdylo Woj06a ( Woj06b . 
The Appendix contains the first few coefficients of the asymptotic expansion under the most general 
hypotheses, then again with the assumption that the minimum is nondegenerate. 

2 The general case. 

In this section, we state and prove our most general result: expressions for the coefficients which appear in 
the asymptotic expansion of integrals of Laplace-type with no smoothness or nondegeneracy assumptions 
on the functions appearing in the integrand. As mentioned in the introduction, the price we pay for this 
generality is that the coefficients are expressed in terms of certain integrals which may be difficult to evaluate 
in practice. In Section[3]below, we discuss various additional hypotheses which are sufficient to alleviate this 
problem. Our main result, Theorem 12.11 below, is essentially a modification of an existence theorem of Fulks 
and Sather |FS61j using a variant, Theorem I2.2[ of a result on inversion of series due to Frame |Fra57j . 

Let x = (x 1 , . . . , x d ) be coordinates on R d . Denote by S^ -1 = {|x| = 1} C K d the unit sphere and 
introduce spherical coordinates p := (x 1 ) 2 + ■ • • + (x d ) 2 and SI = x/p e S d ~ x . Our main result is the 
following theorem. 

Theorem 2.1 Suppose f and g are measurable functions on a measurable set R € M. d which contains as 
an interior point. Suppose further that f attains its unique minimum value of at € R and is otherwise 
bounded away from zero, and that there is a positive integer N and 

1. N + 1 continuous functions fj(il), j = 0, . . . , N with fo(Jl) > such that for some real number v > 

N 

f{p,tt) = p v Y J fiW +o(p N+u ) as p^0, and (2.1) 

3=0 

2. N + 1 functions 9j(0,), j — 0, . . . , N such that for some real number A > 

N 

g(p,Q) = p^-^gjW +o( P N+x - d ) asp^O. (2.2) 

3=0 
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Then if there exists fco > such that 



- kof gd d x 



R 



converges, then there exists an asymptotic expansion 

N 



I e~ k fgd d x = VOfc-°' +A)/ " + °(^ (Ar+A)A 0, k 

jR 3=0 



oo 



(2.3) 



where the coefficients are given by 



c = I r ( 2±a 



/o(!i)- w E«-»(")E 



m / 3+A\ /(r) 



m=0 



/o(«) r 



dfi. 



(2.4) 



where 



ni +712 H h?ij-— m 

nv>l. i— 1 r 



?s i/ie sum of all ordered product^ of r elements of {/i(fi), /^(Q), . . . , //v(f2)} suc/i £/iaf i/ie subscripts 
add to m, and (") := a(a — 1) • • • (a — r + 1)/V!. Empty sums are understood to be 1. 

The first few coefficients of Theorem 12.11 are listed in the appendix. 
Remarks. 

1. If we make the slightly stronger hypothesis that g admits an asymptotic expansion of the form 



3=0 



(2.5) 



for some ft > 0, it is straightforward to make the appropriate modifications to the proof to show that 
there exists an asymptotic expansion 



N 

Jli 3=0 



OO 



(2.6) 



2. If we assume that / / p v and g/p are C at the origin, then by continuity we have that 



p=0 



lim 

p^O 



(p")o G" 



-pf2 



lim 

p-s-0 



p /p(-n) V p To 
P 



p=0 



(similarly for g) which implies fj(-ft) = (— l) J '/ 3 -(fl) and <?j(— O) = (— 1)- J 'gj(f2). If we let 

:= / o( n)-^ e »-«(«) E ( ; 

denote the integrand in (12.41) . then a short computation shows that 

ifc(-fi) - (-1)^,(0). 

It follows then that C2j+i = since it is the integral of the antisymmetric function rjj over the (d — 1) 
sphere. 



4 For example, f& = 6/1/2/3 + 3/f/ 4 + /|. 
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Our proof of Theorem 12 . 1 1 roughly follows that of Fulks and Sathcr [FS61 , although we will use a variant 
of a formula for the inversion of a series due to Frame [Fra57] . There are other, equivalent, formulations 
of series inversion, see for example Kamber [Kam46 ; we choose to use Frame's results because they are 
particularly simple and well suited to our purposes. Because the statement of Frame's result that we need is 
slightly different from the original, and Frame's proof actually holds in slightly more generality than it was 
originally stated, we include the result and the relatively short proof here for completeness. 

Theorem 2.2 JJraJTt Thm. 1] For any real numbers v ^ and q > 0, let u = f(p) and p = f (u) be 

inverse functions defined for p near by the convergent power series 



3=0 



p v y ajp 1 , with oq > 0, and (2-7) 



P q 



'J2 b iU j - (2- 

j=0 



Then the coefficients bj in the inverted power series \2.8\) are given explicitly in terms of the coefficients aj 
of \2. 7p by the inversion formula 

h = -^- q -o u+q)/u t (~T ) a °^ r) forj > °» (2 - 9) 

where empty sums are understood to be 1, (z) ■— a(a — 1) • • • (a — r + 1)/H, and 

(r) _ 

Oj > a ni a n2 ■ • • d nr 

m+«2H h«r=j 

ni>l, i— l.....r 

is the sum of all ordered products of r terms of the set {a\, 02, 03, ... } in which the sum of subscripts is k. 

Proof. First, substituting (|2.7[) into (|2.8p and matching the coefficients of p q shows that 6 = a Q Q ^ ■ ^ n 
particular, since ao > we also have 6 > 0. Choose a circle C in the complex p-plane centered at so that 
the image C under the map p 1— > u(p) is a simple curve which winds around once in the complex it-plane. 

Choose C to be small enough that YljLi a iP' < a o f° r P € C, Yl'jLi b i u ^ 
no other zeroes of ~Y^T=\ a iP l inside C or zeroes of ^j u ^ inside C . 

We can now use Cauchy's integral formula^ to extract the desired coefficients bj\ 



< bo for u € C , and there are 



2m Jqi \u 
If pi du 



2m J c u<i+3+ 1 dp 



dp 



2m J c dp V Q + J J 

= 7T- / — —qp q dp 

2m J c q + j 

where we made a change of variables in the first line, the second line is valid because we assume q > 0, and 
the last line is integration by parts. 

Next, u = p(ao + a(p)) 1 ^ , where at(p) := Y^jLi a jP^ so we have 

h i = h I - J ^p- q - j (ao + a(p)T (q+j)/v p q - 1 dp 
2m Jc q + 3 



1 



^El q } 3 )% {q+j)/ "- r <*(p) r dp 



2m J c q + j f^ Q V r 



00 



q+j ^— ' 



5 Recall that f c z k 1 dz = 2ni if k = and = if k is any nonzero integer. 
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By Cauchy's integral formula, 5^7 Jc aT P ~* 1( ^p = a^p is the sum of all ordered products of r elements of 
the set {cii, 02, . . . } whose subscripts add to j. This sum is only nonzero for 1 < r < j, so we obtain 



The results of Fulks and Sather from which Theorem l2.1l is built are summarized in the following lemm^ 
(for completeness, we include a proof which outlines their argument). Define 

N 
3=0 

and 

R t := {x : f(x) < t}. 

Lemma 2.3 L r FS61] Under the same hypotheses as Theorem \2.1[ for every sufficiently small a > 

e~ kf g d d x = k f e- kt G(t)dt + o{k- {N+x ^ v ) 1 k -> 00 

R JO 



where 



G(t):= / gd' 
'Rt 



Remark. As noted in the remarks following Theorem 12.11 if we make the slightly stronger assumption 
(12.51) . then we can conclude that 

e~ kf g d d x = k f e- kt G{t)dt + (k-^ N+x+mm ^^^ u ), k -> 00 

R JO 

Proof. Fulks and Sather FS61, pl88-189] show that we may assume without loss of generality that g > 0. 
Also, if f R e~ k f g d d x converges for fco > 0, then it converges for all k > ko since e~ k °f g is a dominating 
function for the integrand. We henceforth assume k > fco- 
Given e > 0, define the functions f+(p, f2) and f-(p, f2) by 

N 

f±(p,Sl) :=p v ^2f j W±sp N + v , 

and let 

Rf := {x : /±(x) < t}. 
Choose po > so small that for < p < po 

1 



2. 



/(p,o)-^Ef= /iW 



3. f±(p, £1) is increasing in p for each O, and 

4. i? := {p < Po} C R. 

The first two conditions are obtainable from the assumed asymptotic series for / and g as p — > 0, the 
third condition is insured by the assumption that /o > away from 0, and the last condition follows from 
the assumption that is a point in the interior of R. 



6 Fulks and Sather use this to show that an asymptotic expansion of the form 112.31 1 exists, and they compute the first term 

Co- 
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Then 

I(k) := / e- kf gd d x 
Jr 

can be written as 

/(*)=/ e~ kf gd d x + [ e- kf gd d x. (2.10) 
Jro Jr\Rq 

Since we assume that, away from 0, the function / is bounded away from 0, there exists A > such that 
/ > A for p > po, which implies the second integral above is 0(e~ kA ) as k — > oo. 
Next, define 



and 



I±(k) := / e- k f±gd d x, 
JRo 

G±(t) := / gd d x. 

JR± 



Since /_ </</+, we have 

RfdRtdRi- (2-11) 

Choose a > so small that i?~ C Rq. Combining jFS61[ Lemma 2] with the same argument we used above 
to conclude that the second integral in (|2.10[) is exponentially small, we obtain 



I±(k) = [ e- kf± gd d x+ [ e- kf± gd d x 
Jr£ JRo\Ri 

= [ e- kt dG ± (t) + 0(e- kA± ) (2.12) 
J o 

= k I e- kt G±(t)dt + 0(e- kA± ) (integration by parts) 
Jo 

for some positive constants A±. 
By (|2.1ip we have 

pa pa pa 

e- kt G+(t)dt< / e - kt G(t)dt< / e- kt G-{t)dt. (2.13) 



We also have (by the definition of f±) 

I+(k) < f e- kf gd d x < I_(k) 
JRo 

which, by (|2.12[) . and (|2 . 10[) and the subsequent argument, gives 

e- kt G+(t)dt + 0(e- feA +) < I(k) + 0{e~ kA ) <k P e - kt G-(t)dt + 0(e~ kA -). 

Jo 

Subtracting (|2.13[) then yields 

/•a 

e- kt {G+(t)-G-(t))dt + 0(e~ kA +)<I(k)-k e - kt G(t)dt + 0{ e - kA ) (2.14) 

Jo 

< k I" e~ kt {G^{t)-G + (t))dt + 0{e- kA -). 
Jo 

Fulks and Sather |FS61[ pl91] show that there is a constant c > such that 

k P e- kt (G^(t) - G+(t))dt = eck-^+^l" + o(k-( N+x V»). 



'In fact, c ■ 



I r (^) fs*-i 3o(n)/ (n)-i-(^+A)/-da. 
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Equation (|2.14[) thus becomes 

I(k) - k 

Hence we conclude that for every e > 



-kt 



G(t)dt 



lim 



i(k) 



-kt 



G{t)dt 



k (N+\)/u < 



sc. 



from which the lemma follows. 



Our proof of Theorem 1 2. II is now essentially an estimate of the function G(£) using the inversion formula, 
Theorem 12. 2[ of Frame. 



Proof of Theorem 12.11 As in the proof of Lemma POl we henceforth assume that g is positive and that 
k > ko so J R e~ k f g d d x converges. By Lemma 12.31 it is enough to show that for a > sufficiently small, 
there exists an asymptotic expansion 



-kt 



G(t)dt 



N 

££fc-C*+A)/* + (fc-("+A)/»') j k 
3=0 



oo 



where the coefficients Q are given by (I2.4j) . To this end, choose po small enough that / is increasing in p for 
< p < po, and then choose a small enough that R a := {x : /(x) < a} C {p < po}. Then for each £ with 
< t < a, the equation t = f(p, has a unique solution p(t, f2) for p which is continuous in fi. Substituting 
the series (12 .2|) for g and performing the integration over p yields 



G(t) 



s-*- 1 Jo 



p(t,n) 



g{p,n) P d - 1 dpdn 



S d~i 



N 



E#tp j+a (*^)+ o (^ +a ) 



(2.15) 



Now we use Theorem l2.2l to explicitly estimate the powers of the inverse function p(t, thus obtaining 
explicit estimates for G(t). 
Let 

_ (fj(Sl) for < j < N, and 
flj '~ [0 forj > N. 

Then t:=f = p» E^o^'- 

Let us briefly consider the inverse function of £ = /(p, O), for fixed O, in more detail. Let u = t 1 /" , so that 

/o(f2) 1 /", which is positive by assumption. 



l/u 



= p \ ^27Lq &ifP ) . One may compute that j- 



dp 



p=0 



Moreover, u is an analytic function of p in a neighborhood of 0. Since u(p = 0) = 0, the inverse function 
is also analytic in a neighborhood of 0, and we may express it as a series p = u^™ bjU J . Taking positive 

powers of this series yields series p q = u q 2j*lo b^u 3 . In particular, the two series u u — p v Xj*lo a iP > an< ^ 

p q = u q J^'jLo tfj " J satisfy the hypotheses of Theorem 12.21 

Hence, with u = i 1 /", Theorem 12.21 yields p q = u q J2j^-ob^u J , wnere the coefficients are given 

by (|2.9I) . Clearly, depends only on {do, Oi, . . . , dj}. Moreover, the remainder after N terms, 

u q J2N+i b^ityu^ is uniformly bounded for Q € S^ 1 and < u < a 1 / 1 '; that is, 



Since u — t 1 '", we get 



j=N+l 



bfiny = o(u 



N 



N+l+q\ _ 



) = o(u N+q ). 



p q = ^"^^ ) ( fl ) ti/ "+ ( f(W+,)/ ")' 
J=0 



(2.16) 
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Substituting (12.161) into the expression (|2.15l) for G(t), we obtain 

N N 



j=0 1=0 

After some rearrangement H one obtains 



1 



S d-i 



N + A 



N 



where 



3=0 



V3 = 12j^j 9 3-mr i+x) m d n. 

1=0 J b 



Finally, recall that for s > — 1 



k e - kt t s dt = k- s r(s + i) + o(k-°°). 

Jo 



We multiply G(t) by e and integrate termwise to get 



k [ e- kt G{t)dt = k ( e~ kt (V^t^+o^) | dt + 0(/c~°°) 

7o 70 v=° / 



A' 



5^( 3 .jfe-a+W" + o(A-( w + A )/") 



where, since T(s + 1) = sr(s), 



m=0 r=l V ' Jl)l ' 



<K1 



(2.17) 



(2.18) 



We conclude this section with a short discussion of the combinatorial quantities which arise in Theorem 
2.11 The quantities 

fm) = y , fnxfni ' ' ' fn r (2-19) 



ni+ri2^ \-n r =m 

rii>l, i=l ,r 



also arise in the 1-dimensional case studied by Wojdylo in |Woj06a , |Woj06bj, where he calls them "partial 
ordinary Bell polynomials." They are relatively simple combinatorial objects to compute; for example, f„^ 
may be computed as the coefficient of x n appearing in (fix + fix 2 + + • • - ) r . The following result 
describes a simple recursive algorithm for their computation. Explicit algorithms suitable for computer 
implementation may be found in |Woj06a| , |Woj06b| . 

Proposition 2.4 For n > we have fn — fn> an d f or I < r < n, the coefficients /„ can be computed 



recursive 



ly via /« = J2 U-ifj 

2=r— 1 



(r-1) 



8 To obtain (J27THJ, define hjj := / 5<J _i g 3 (n)b\ j+X) {ft)dfl. Then the double sum in iTTTTt is J2jLo T.1L0 ' l J>i* i_ ^ _ • ° f 
course, if both j and I are large, then the corresponding term is absorbed into the error term. So in fact we only need consider 

j = 0I+j<iV 3=0 i = m=0 i+j=m m=0 (=0 

Putting h m _i : i into this expression then yields 112.181 1. 
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3 Further Simplifications 



In this section, we make additional hypotheses which allow one to evaluate the integrals which appear in 
Theorem 12.11 In Section I3TT1 we assume first that /o(^) is constant, that is, that f/p 2 is continuous at the 
origin, which allows us to make slight simplifications to the coefficients £j . We then assume that / and g are 
smooth enough to admit at least finite order Taylor series expansions, and evaluate the integrals appearing 
in the coefficients £j of Theorem 12.11 explicitly, thus yielding an asymptotic expansion of J R e~ k ' g ei d x in 
terms of the Taylor coefficients. In Section 13.21 we make the assumption that / attains a nondegenerate 
minimum at (this is what is traditionally assumed to make Laplace's approximation) and also that / and 
g are smooth enough to admit at least finite order Taylor series expansions. These hypotheses also allow us 
to evaluate the integrals appearing in the coefficients £j explicitly, and lead us to a proof of Theorem 11.11 
which was stated in the introduction. 

3.1 Simplifications with /o(fi) constant, and with enough smoothness for Taylor 
series. 

First, we simplify Theorem 12.11 under the assumption that fa(0.) is constant. Then, we show that if fa(£l) is 
constant and / and g are given as Taylor series with remainder, the integrals appearing in the coefficients 
Q in Theorem 12.11 can be explicitly computed to obtain a closed formula for Q involving the coefficients of 
the Taylor series. 

If fo(£l) = fo is constant, then the coefficients of Theorem l2 . 1 I simplifv significantly since the denominators 
consist entirely of terms involving fo(Cl). In general, one does not expect fo(£V) to be constant, but in the 
common case that v = 2 (discussed in l3.2l below). an appropriate coordinate change guarantees fo{&) = 1/2. 

Before stating the main asymptotic expansion with the additional hypothesis that fo is constant, we 
mention an easy lemma characterizing the condition that fo is constant. 

Lemma 3.1 Suppose f is a measurable function on a measurable set R C M. d containing as an interior 
point, such that for some positive integer N and some positive real number v 

N 

f( P ,n)= P »Y,fiW +o( P N +n. 

3=0 

Then p~ v f{p, f2) is continuous at p — if and only if fo(il) is constant. 

Proof. The function p~ v 'f(p,Q) is continuous at p = if and only if c := lim^o p~ P /(p, fi) = /o(fi) is 
constant. ■ 

Corollary 3.2 //, in addition to the hypotheses of Theorem \2.1[ p~ v f(p,Q) is continuous at p — 0, so that 
/o(^) = fo is constant, then there exists an asymptotic expansion 

N 

/ e- fe V d x = y2Qk- (j+X)/ " + o(fc- (W+A)/,/ ), k -> oo (3.1) 

jR 3=0 

where the coefficients are given by 

1 / \ 3 m / j+X \ r 

Ci = ~r (J±*) f^ +x),v EE ; )fo r / 93-m(n)f^(n)dn. 

m =0r=l ^ ' JSd 1 

Empty sums are understood to be I. 

Now we suppose further that / and g are smooth enough to be given as Taylor series (with remainder), 
and give closed form expressions for the coefficients appearing in Theorem 12.11 (still with the assumption 
that /o(O) is constant); that is, we explicitly evaluate the spherical integrals appearing in the coefficients Cj 
of Theorem 12. II to obtain formulas in terms of the Taylor coefficients. 



11 



To relate the Taylor series expansions with the radial expansions of Theorem 1 2. 1[ we introduce spherical 
coordinates (p, 0i, 02, ■ ■ ■ , 0d-i) on R d ; these are related to the Cartesian coordinates x = (a; 1 , . . . , x d ) via 



x 1 = p cos 0i, 

x 2 = p sin 0i cos 02 , 



x d 1 = p sin 0i sin 02 • • • sin 0^-2 cos 4>d-i, 
x d = p sin 0i sin 02 • • • sin 0^-2 sin0^_i. 



(3.2) 



The coordinate ranges are < p < oo, < 4>j < tt for j — 1, . . . , d — 2 and < (j>d-i < 2ir. The solid angle 
element in spherical coordinates is 

dQ = sin d ~ 2 0i sin d ~ 3 02 ■ ■ ■ sin0d_2d0i ■ ■ ■ d'Pd-i- (3-3) 

According to (J3T2J) , the unit vector f2 = x/p may be expressed as 

= (cos 0i , sin 0i cos 02 , . . . , sin 0i sin 02 • • • sin <fid-i) ■ 

Suppose / and g are given as Taylor series with remainder, in standard multi-index notation, as 

N+v 

/(x) = V -D a f(0)x a + R N+v+l (x), (3.4) 

\a\ — v 

for some integer v > (actually, if / is to have a unique minimum of at 0, the parameter v must be an 
even integer > 2), and 

N 1 

■9« = E o^s(0)x" + J4+i(x). (3.5) 

1/9 1=0 P ' 

Observe that (|3.5[) implies A = ef. To compare these expressions with the radial power series used in Theorem 
12.11 we write (|3 .4[) and (|3.5[) in spherical coordinates to obtain 

/(p,Q)=Xy E ^ffi^ + 0(p^ +1 )and 

j=i/ \ a \=3 
N 



^ ) = E^E^r flo + °^ +I ) 

from which we see that 



a! 

J=0 \ a \=j 



m= E ^r 1 ^-d ft (o)=5:^)o« (3.6) 



Corollary 3.3 //, m addition to the hypotheses of Theorem \2.1\ the functions f and g are smooth enough 
to be given as the Taylor series and \S.5\) where fo(£l) — fo is constant, then there exists an asymptotic 

expansion 



3=0 

where the coefficients are given by 



N 

/ e-^g d d x = k~ d ' v V Qk- j/v + 0(k^ N+d+l ^ 1 ') 

JR 



rn—0 r—1 ^ ' 



x E E E ••• E vp+«r+~+*r R\ ai \... a ! 

\fi\=j-7n (n 1 -\-7i 2 -\ \-n r =m\ \ax\=m+v |a r |=n r +^ 

\ ?ii>l, i—l,...,r j 
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where for a multi-index a — (a 1 , a 2 , . . . , a d ), 



\dli f I if d even 1 (a — 1)!! . . 

y Q := (2vr d/2 <^ rr-f- ., , } -r ,, even a . 

v ; \ V2/tt i/d odd J (|a| + rf - 2)!! v ; 



(3.7) 



Empty sums are understood to be 1. 



Recall that for an integer n, the double factorial is defined^! to be nil := n(n— 2)(n — 4) • • • (n— 2\n/2\ +2) 
where [n/2\ is the greatest integer less than n/2, that is, 



n(n — 2) • • • 4 • 2, if n even 
n(n — 2) • • • 5 • 3, if n odd 



and we adopt the convention that (—1)!! = 1. Moreover, for a multi-index a, we dehne all := a 1 !!a 2 !! • • • a d \\. 
Finally, even(a) := 1 if all components of a are even, and otherwise. 

Proof. By Corollary 13.21 and Remark 1 following Theorem 12.11 we need only compute the integrals 
Joi-i 9j-mfm dQ. The coefficients /j(f2) and gj(Cl) obtained from the Taylor series (J3T6J) and p.5j) are 
given by (|3.6|) . From (|2.19[) , the integrand of interest is therefore given by 

»_(n)jS>(n) = n; n; e - E 

|/3|=j— m f ET=l Ti * =m "\\ai\=ni+v \a r \=n r +v 
\rii>l, i— l,...,rj 

Thus, it remains to compute the integral 



j3\a x \---a r \ 



(3.. 



S d-i 



for an arbitrary multi- index a. 

Let a = (a 1 , ... , a rf ). By (072]) and ([33]) . we have 



,or+crH ha a +d-2 



n^dfi = ( cos Q tpi sin 1 * 1 " 1 ~ <^J X ^cos" </?2 sin' 

: (cos"" 1 <^ rf _i sin Qt! ^-lj ^i • • • dtpd-i 



a 3 +a i A \-a d +d-3^ 



(3.9) 
(3.10) 

For nonnegative integers a and b, an elementary computation (using reduction formulas, for example [GR961 
2.510]) shows that 

6 , (o-l)!!(6-l)!l , 
sin a; cos xdx= ; ttt; c\(a,b), 



where 



and 



where 



{a + b)\\ 

if & is odd, 
ci(a, b) := ^ 7r if both a, b are even, , 
2 if a is odd and & is even 



2n ■ a b . (o-l)!!(6-l)H , 
sin xcos xax — ; — — C2(a,b), 



c 2 (a,6) 



(a + 6)!! 

if a or (or both) is odd, 
2tt if both a, & are even, 



9 We could define w a by the nonzero formula appearing in 1 13. 71 1 if we adopted the nonstandard convention that for n even, 

nil = n(n - 2) • ■ ■ 4 • 2 • = 0. 

We will not do this. 
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Combining this with (|3.9p shows that 



d-i (a 1 - 1)!! {Y? =1+x a p + d-2- 
n a dn=]J N — x II c i E o? + d-l~l,a l \ xc 2 (a d . u a d ). 



qd-1 

'=1 



E^^ + d-l-nH 




Now, if any a p is odd, the entire integral is zero, so we assume for the moment that all a p are even, which 
implies aP is even. If d is even, then ci ^X)p=;+i a p + d — I — I. a l ^j = tt if I is odd and 2 if Z is even. If 

d is odd, then the opposite is true, that is, c\ ^X)p=;+i ck p + d — I — I, a l ^j = 2 if I is odd and n if £ is even. 
In any case, 02(0:4—1,0:4) — 2n. It is now a simple matter to count the factors of 2 and tt appearing in the 
above product in each case of d even or odd, which yields the theorem. ■ 

3.2 Assuming a nondegenerate minimum. 

In this section we consider the case that / has a nondegenerate minimum at 0. It follows (from the proof 
of Theorem ll.il below) that v — 2 and that fa(SY) > 0; in fact, by an appropriate choice of coordinates, we 
can guarantee that fo(0.) = 1/2. Hence, Corollary 13.21 applies and we obtain Theorem ll.il which was given 
in the introduction. 

Proof of Theorem 11.11 Since / has a nondegenerate minimum at 0, the Hessian can be diagonalized, 
Hf(0) = Q~ l DQ, for some orthogonal matrix Q and some diagonal matrix D. The eigenvalues of Hf(0) 
are all positive, so we can unambiguously define \pD. Let P := y/DQ, and define y := Px. Denote the 
Hessian of / with respect to the new coordinates y by H y f. Then H y f(0) is the identity matrix. Moreover, 

the Jacobian of the coordinate transformation, since Q is orthogonal, is det VD = J H x f(0) . Thus, the 
integral of interest becomes 

/ e- fe ^.g(x)d d x= y/detHsHO)' 1 f e~ k f^g(y)d d y. 

JR J PR 

Now we apply our results to the last integral above; in particular, we introduce spherical coordinates 
p = |y| and Q — y / p. Again since / has a nondegenerate minimum at 0, the first term in the radial expansion 
of / at p = is 

\a\=2 3=1 

whence we can apply Theorem 13.21 On the other hand, the smoothness assumptions on / and g imply that 
we can express them as Taylor series with remainder as 



/(y) = E ^ D yf(°)y a + °(p N+3 ) and 9(y) = E ^ D t^°)y a + o( P N+1 ), 

\a\=2 ' \a\<N 

and we can apply Corollary 13. 31 to obtain 

Q = ^/det^/CofV to) 2^ V V f-^V 



m— r— 1 



r I 



^ ^ \ ... \ Wi3 +ai + ... +ai 



DPg(0)D?f(0)---D?f(0) 



^ ^ ^ ^ -P+ax+-+«r p\ ai \...a r \ 

\P\=j—m fn 1 +n 2 -\ \-n r =m\ \ ai \= ni +2 \a r \=n r +2 

\ ni>l. i— l,...,r / 

Now, if j is odd, then \/3 + oti + • • • + a r \ = j — m + m + 2r = j + 2r is odd, whence wp+ ai -\ |_ Qr = since 

at least one component of (3 + ot\ + ■ ■ ■ + a r must be odd. This means C2.7+1 = and we can replace j by 2j. 
Inserting the value of wp+ ai ^ y ar and using the fact that for each positive integer j we have 



r (| + j) = 2-^-V- d ' 2 (d + 2j - 2)!! j 



1 if d is even 
/tt/2 if d is odd 
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the coefficients ^j, after some simple cancelations, become 



>2j 



(2^/2 ^^ orM | +iA (d + 2i-2)!! 



EE 2 ' 



Vdetflx/CO)^,^ V r J (d+2j + 2r-2)U 



x E E E ••• E 

\f)\=2j-rn (ni+Tl2-\ \-n r =m\ | Ql |=rii+2 \a r \=n r +2 

\ ni>l, i— l,...,r / 

x (P + «! + ■■• + a r )U even(/3 + «! + ••■ + a,,). 

Finally, one may compute that 

(d+2j-2)\\ _(-iy 



2 

to obtain the theorem. ■ 



r 



(d + 2j + 2r-2)\\ r! 



The coefficients Q of Theorem 11.11 can be written in terms of the derivatives of / and g with respect to 
the original variables x by inserting the appropriate entries of P _1 . 
Finally, we relate Theorem 1 1.1 1 to Laplace's approximation. 



Corollary 3.4 (Laplace's Approximation) Let R C M. d be a measurable set which contains as an interior 
point. Let f be a twice differentiable function on R which attains a nondegenerate, unique minimum value 
ofO at 0, and let g be a continuous function on R. Assume moreover that there exists a positive real number 
ko such that 

e~ kof g d d x 



converges. Then 

[ e-Vg rf d x = ( ^ - g(0) + 0{k- d '^), k -+ oo. 

Proof. Since g is continuous, we have A = d and go(£l) = g(0) is constant. The Corollary then follows from 
the formula for Co of Theorem 11.11 ■ 



3.3 The 1-dimensional case. 

We compute the coefficients of Theorem 12.11 in the 1-dimensional case, where several simplifications occur, 
thus rederiving the recent results of Wojdylo Woj06a , |Woj06b] Thm 1.1]. 

Corollary 3.5 Suppose that for x — > 0, we have f(x) = x v X^jLo a o x ^ °{x N+v ) for some v > and 

9( x ) = S^o^i 2 '"' +°{x N ). Then if there exists ko > such that f_ a e~ k °J g dx converges, there exists an 
asymptotic expansion 



f e'^gdx^k- 1 ^ £ C 2j -fe- 2i/v 
where the coefficients are given by 



o(fc -(2LJV/2j+l)/ V j j k _^ QQ 



( 2j m , 2j + l \ ^ 



in which am is the sum of all ordered products of r terms of the set {a\, a 2 , ... } such that the subscripts add 
to m, \_N/2\ denotes the largest integer less than N/2 and empty sums are understood to be 1. 
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Proof. Note that the hypothesis on g implies A = d in Theorem 12. II 

Introduce "spherical coordinates" p(x) — \x\ and Q(x) = xj \x\ = ±1 on R. Note that 5° = {±1}- Since 
= 1 if j is even and = fl if j is odd, 

/ tt j dn = 2ej (3.11) 
Js° 

where £j = 1 if j is even and otherwise. 

The series for / in polar coordinates is f(p,Q) — p v Y2jLo p* ■, so that fj(il) — djiV , which is equal 
to a,j if j is even and a^fi if j is odd. Similarly, <?j(r2) = bj if j is even and bjVL if j is odd. Using these facts, 
and also the observation that the power of 51 in fm is m, we obtain from Theorem 12.11 that 



' r 

r—1 



a a m 



from which it is clear that Q — for j odd, whence we obtain the desired result. ■ 

It is worth mentioning that in |Woj06a| Sec. 6.1], to check the correctness of his formulas (essentially the 
same as Corollary [33]), Wojdylo applies his expansion to the T-function by writing 

/>oo />oo 

k\ = T(k + 1) = J e-h k dt = k k+1 e- k J e' k{x " ln{x+1)) dx. (3.12) 

The resulting expansion (i.e., applying Corollary 13.51 to (|3.12|l ) can be used to obtain the following beautiful 
expression for the Stirling series (cf. [AS641 6.1.37], [Sti30j ) 

ik -k /tt/, 1 1 139 571 163879 \ , 

fc ~ k k e k v2wk 1 H 1 ^ r 7 H F H , k -> oo. 

V 12k 288fc 2 51840fc 3 2488320fc 4 209018880/c 5 / 



Corollary 3.6 TTie Stirling series is given by 

2j 



j=0 r=0 



k\ ~ k k e~ k V2^k S ^~ (2j + 2r ~ a 2? 



where at' is i/ie sum o/ all ordered products of r terms of the set {ai, 02, ... } in which the subscripts add to 
m, where a,j = j^, and where (—1)!! is understood to be 0. 
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4 Appendix 

Listed below are the first few coefficients appearing in the asymptotic expansion, Theorem 12.1 



V \ V 



gi(») (A+D A(n)go(n) 



+i 



, =lr ,Mj g2(») A+2 g0(")/2(")+gi(")/l(") | f-^\ 90(»)/l(O) 2 , n 



-r , 

-^\ gi(n)/i(n) 2 + 2 g0 (n)/ 1 (n)/ 2 (n) f-^ ^Wjf .n 

v 2 ; / (o)^+ 2 i 3 ; /o( ^+3 > 

L T (h±±) [ g 4 (») A+4 go(n)/ 4 (n) + gi(n)Mn) + g2 (n)/ 2 (») + g3 (n)/i 
-A±i\ /a( n)2 + 2 g0 (n)/ 1 (n)/ 3 (n) + 2 gi (»)/ 1 (»)/ 2 (o) + g 2 (»)/i(») 2 

2 / /o(^+ 2 

^ 3, 9 o(»)/ 1 (») 2 / 2 (»)+,g 1 (»)/ 1 (0) 3 ^ f- .g (»)A(») 4 ^ 



3 / /o(fi)^ +3 V 4 y /o (f7)^+4 

Next, we give the first two nonzero coefficients appearing in Theorem 11.11 Here, the coefficients are 
expressed in terms of the derivatives of / and g at with respect to y = Px, where P is the matrix such 
that H y f(0) = Id. In particular, if H x f(0) = Q~ X DQ for a diagonal matrix D and orthogonal matrix Q, 
then P = \J~DQ. Here, Q = for odd values of j (see Remark 2 following Theorem 12.11 or the proof of 
Theorem 1 1.1|) . and we give only the nonzero values. 

(27r) d / 2 o(0) 

Co = — (Laplace's Approximation) 

y/detH x f(0) 

— ^-(a - l)Heven(a) + ~ — - — ; — v - — (a\ + a 2 - 1)!! even(ai + a 2 ) 

£ — ' a! z — ' ai!a 2 ! 

|a|=4 |ai| = |a 2 |=3 / y=0 

Finally, we give the above coefficients precisely for small d. (Recall that for d = 1, our results coincide with 
Corollary [331 originally due to Wojdylo |Woj06a , Woj06b].) 
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Dimension d = 1 



Co = 

(2 = 



2; f \g"{Wm 2 V(°)/ (3) (o)/"(o) - I g( o)/( 4 )(o)/"(o) + A. g( o)(/( 3 )(o)) 2 



C 4 = (/ ,,^ 13/2 (^ (4) (0)(/"(0)) 4 - ^.9 (3) (0)/^(0)(/"(0)) 3 

+ ^"(0) (-^/ (4) (0)(/"(0)) 3 + ^(/(3)( )) 2 (/"(0)) 4 ) +g'(0) r-I/W( )(/"(0)) 3 + |/W( )/( 3 )(0)(/"(0)) 2 - |(/( 3 )(0)) 3 .r(0) 



+ 



fl(0)(- ^/ (6) (0)(/"(0)) 3 + ^(/( 4 )(0)) 2 (r(0)) 2 + l/(5)( )/( 3 )(0)(r(0)) 2 - |/W(0)(/( 3 )(0)) 2 /"(0) + ^(/ (3) (o))' 



Dimension d = 2: 

2tt 5 (0) 



Co- 



^dct #*/(()) 



^ + d^g)) - \d v ,g (d yl d 2 y2 f + d 3 yl f) - \d y ,g (d^f + d y2 f) 



s/AetH x f(0) 



Dimension d = 3: 

Co = 

c 2 = 



(27T) 3 / 2 g(0) 
y/dctH x f(0) 

(2tt) 3 / 2 
^/dct #*/(()) 



\ {dyig + d 2 y2 g + d 2 y3 g) - \d y ig (<9 3 i/ + dyid^f + d y id 2 y3 f) - \d y2 g (d^f + d y2 d 2 y3 f + d 2 y id y2 f) 



- \dysg {d 3 y3 f + d 2 y2 d y3 f + d 2 yl d y3 f) 

+ I {(9 2 yl d y2 f) 2 + (d 2 yl d y3 f) 2 + (d 2 y2 d y3 f) 2 + {dy^ff + (d y ,d 2 y3 f) 2 + (d y2 d 2 y3 .f) 2 ) 

+ \ (dy^ + dy^f) + d^f (d 2 yl dy 2 f + dy.d 2 y3 f) + d 2 y3 f (d^f + d 2 y2 dy 3 f)) 

+ \ (<V%/ d v ld yi + d l* d v 3 f d yi d v 3 f + d y ^d 2 3 f d^d^f) + d y i d y 2 d y3 /) 



y=o 
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